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POINTWISE BEST APPROXIMATION RESULTS FOR GALERKIN FINITE ELEMENT SOLUTIONS 

OF PARABOLIC PROBLEMS 

DMITRIY LEYKEKHMANt AND BORIS VEXLER* 

Abstract. In this paper we establish a best approximation property of fully discrete Galerkin finite element solutions of second order 
parabolic problems on convex polygonal and polyhedral domains in the L°° norm. The discretization method uses of continuous Lagrange 
finite elements in space and discontinuous Galerkin methods in time of an arbitrary order. The method of proof differs from the established 
fully discrete error estimate techniques and for the first time allows to obtain such results in three space dimensions. It uses elliptic results, 
discrete resolvent estimates in weighted norms, and the discrete maximal parabolic regularity for discontinuous Galerkin methods established 
by the authors in nu. In addition, the proof does not require any relationship between spatial mesh sizes and time steps. We also establish an 
interior best approximation property that shows a more local behavior of the error at a given point. 

Key words, parabolic problems, finite elements, discontinuous Galerkin, a priori error estimates, pointwise error estimates 
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1. Introduction. Let n be a convex polygonal/polyhedral domains in R N , N = 2,3 and / = (0, T). We 
consider the second order parabolic problem 

d t u(t, x ) — A u(t, x ) = f(t, x), (t, x) £ I x fl, 

u(t,x) = 0, ( t,x)£lxdfl , (1.1) 

it(0, x) = uo{x), ie!l. 

For the purpose of this paper we assume that / and uq are such that the unique solution u of (11.Il l fulfills u £ 
C(I x G) H C(I ; Hq (fl)). To achieve this, we can for example assume that the right-hand side f £ L r (I x Q) 
with r > ^ + 1 and uq £ C(S2) fl Hq (G), cf., e. g., 1 421 Lemma 7.12], but other assumptions are possible. 

To discretize the problem we use continuous Lagrange finite elements in space and discontinuous Galerkin 
methods in time. The precise description of the method is given in Section [2] Our main goal in this paper is 
to establish global and interior space-time pointwise best approximation type results for the fully discrete error, 
namely, 

T 

\\u-u k h\\ L °°(ixn) < C|ln/i|ln-||u-xllL“(/xn), (1-2) 

where Ukh denotes the fully discrete solution and x is an arbitrary element of the finite dimensional space, h is 
the spatial mesh parameter and k stands for the maximal time step. Such results have only natural assumptions 
on the problem data and are desirable in many applications, for example in optimal control problems governed by 
parabolic equations. 

Most of the work on pointwise error estimates for parabolic problems were devoted to establishing optimal 
convergence rates for the error between the exact solution u{t) and the semidiscrete solution ?//, ( t ) that is contin¬ 
uous in time, 00 0 0 [201 E2 [231 EH HOI [32] [33] EH . The best approximation results for the semidiscrete error 
u(t) — Uh(t ) in L°°(I x G) norm can be found, for example, in Ifl4ll32l . 

Results on fully discrete pointwise error estimates are much less abundant. Currently, there are several tech¬ 
niques available for obtaining fully discrete error estimates. One popular technique splits the fully discrete error 
into two parts as u — Ukh = (u — Uh) + [uh — Ukh )■ The first part of the error is estimated by the semidiscrete error 
estimates and the second part of the error is treated by using results from rational approximation of analytic semi¬ 
groups in Banach spaces. Thus, for example, optimal convergence rates for backward Euler and Crank-Nicolson 
methods were obtained in |]33l (see also 1 401 Sec. 9] for treatment of general Pade schemes). A similar technique 
uses a different splitting, u — Ukh = (u — RhU ) + (RhU — Ukh ), where Rh is the Ritz projection. In this approach 
the first part of the error is treated by elliptic results and the second part of the error satisfies a certain parabolic 
equation with the right-hand side involving (it — ///,'«), which again can be treated by results from rational ap¬ 
proximation of analytic semigroups in Banach spaces ED (see also ll40l Thm. 8.6]). For smooth solutions, both 
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approaches above produce error estimates with optimal convergence rates. However, in many applications these 
two techniques require unreasonable assumptions on the data, as well as on the regularity of the solution. As a 
result, the best approximation property (11.2b can not be derived, except for the one-dimensional case 1431 . 

Another approach, that is more direct, is based on the weighted technique. For N = 2 and low order time 
schemes, this technique works rather well and allows one to obtain sharp results. Thus, in j£[ (see also (25] Thm. 
4.1]) optimal convergence error estimates of the form 

\\u(t n ) -u kh (tn) ||z/~(n) < C\lnh\ x™< n ( kq \\ d t u h°°((o,t m )xn) + ^ 2 |l-D 2 w|U»((o,t m )xa)), 

for piecewise constant and piecewise linear time discretizations, i.e. q = 1 and q = 2 , correspondingly, were 
derived on convex polygonal domains (the result in J9] actually holds even on mildly graded meshes). The best 
approximation property of the form ( 11 . 2 I > was derived in ll28l on convex polygonal domains without any unnatural 
smoothness requirements. However, for N = 3, the weighted technique is much more cumbersome and as of 
today, there is no three dimensional pointwise best approximation results or optimal error estimates even for 
backward Euler method. 

In this paper for the time discretization we consider discontinuous Galerkin (dG) methods of an arbitrary 
order. These methods were introduced to parabolic problems in lfl 2 l and deeply analyzed in fTTl . There are a 
number of important properties that make dG schemes attractive for temporal discretization of parabolic equa¬ 
tions. For example, such schemes allow for a priori error estimates of optimal order with respect to discretization 
parameters, such as the size of time steps, as well as with respect to the regularity requirements for the solution 
ME- Different systematic approaches for a posteriori error estimation and adaptivity developed for finite element 
discretizations can be adapted for dG temporal discretization of parabolic equations, see, e. g., 03 EH- Since the 
trial space allows for discontinuities at the time nodes, the use of different spatial discretizations for each time step 
can be directly incorporated into the discrete formulation, see, e. g., GZ). Compared to the continuous Galerkin 
methods, dG schemes are not only A-stable but also strongly A-stable US. An efficient and easy to implement 
approach that avoids complex coefficients, which arise in the equations obtained by a direct decoupling for high 
order dG schemes, was developed in |]29l . 

Our approach in establishing (11.2b for dG methods is more in the spirit of the work of Palencia l26l and does 
not require semidiscrete error estimates or even any error splitting. Moreover, it does not require any relationship 
between the spatial mesh size h and the maximal time step k, which is essential for problems on graded meshes. 

Our approach is based on two main tools: The newly established discrete maximal parabolic regularity re¬ 
sults EGD for discontinuous Galerkin time schemes and discrete resolvent estimates of the following form: 

||(2 + A/ ( )“ 1 x||z,‘~(n) < pj||xlU~(n): forz<G<C\£ 7 , for all % e V* = V h + iV h , (1.3) 

where 14 is the space of continuous Lagrange finite elements and 

£ 7 = { z e C | |arg(z)| < 7 } , (1.4) 

for some 7 £ (0, and the constant C that may contain |ln/i| but must be independent of h otherwise. Such a 
discrete resolvent estimate can be shown directly or by showing stability and smoothing results of the 

semidiscrete solution operator Eh(t) = e~ Aht I20l l32l . The first approach is preferable since it establishes (11.3b 
for an arbitrary 7 £ (0, f), while the second approach via theorem of Hille (see, e.g., Pazy [271, Thm. 2.5.2) only 
guarantees existence of some 7 £ ( 0 , f). 

In this paper we also establish a local version of the best approximation result l ira . This result (cf. Theorem 
12.2b shows more local behavior of the error at a fixed point. For elliptic problems such estimates are well known 
(cf. Il34l 13611441 ). but for parabolic problems the only result we are aware of is in 11281 . which is stated for convex 
polygonal domains without a proof and fl5] [TBj that are global in time. To obtain this result, in addition to the 
stability of the Ritz projection in L°°(fl) norm and the resolvent estimate (11.3b . we need the following weighted 
resolvent estimate 

\\a%(z + At t )~ 1 x||i / 2 (n) < lk^xlU 2 (n), forzeC\S 7 , for all x G (1.5) 

with a(x) = ^J\x — Xq \ 2 + K 2 h 2 . This estimate is established in Theorem 14.II The estimate (11.5b is somewhat 
stronger than the corresponding resolvent estimate in L°° norm, meaning that (11.3b follows rather easily from 
(O (modulo logarithmic term |ln h |), but not vice versa. 
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The rest of the paper is organized as follows. In the next section we describe the discretization method and 
state our main results. In Section |3l we review some essential elliptic results in weighted norms. Section [4] is 
devoted to establishing resolvent estimate in weighted norms. In Section[5J we review some results from discrete 
maximal parabolic regularity. Finally, in Sections[6]and[7] we give proofs of global and interior best approximation 
properties of the fully discrete solution. 

2. Discretization and statement of main results. To introduce the time discontinuous Galerkin discretiza¬ 
tion for the problem, we partition the interval (0, T] into subintervals I m = (< m _ i, t m \ of length k m = t m — t m _ i, 
where 0 = to < ti < • • ■ < <m-i < £m = T. The maximal and minimal time steps are denoted by 
k = max m k m and fc m i n = min m k m , respectively. We impose the following conditions on the time mesh 
(as in ED or [21): 

(i) There are constants c,/3> 0 independent of k such that 

^min ^ Ck^. 

(ii) There is a constant k > 0 independent of k such that for all to = 1,2,..., M — 1 

— 1 ^ km ^ 

K < -- < K. 

Km -\-1 

(iii) It holds k < \T. 

The semidiscrete space of piecewise polynomial functions in time is defined by 

X q k = {u k €L 2 (I;H^n)) \u k \ Im £V q (H^(Q)), m = 1,2,..., M } , 

where V q (V) is the space of polynomial functions of degree q in time with values in a Banach space V. We will 
employ the following notation for functions in Xj* 

4 = lim u(t m + e), u“ = lim u(t m ~e), [u\ m = - u“. (2.1) 

£—^0 ' £—^0 ' 

Next we define the following bilinear form 

M M 

B{u,ip ) = '^2{d t u,(p) ImX n + (Vu, V ^)/ Xf 2 + Y (Mm-i, + (“o', ¥?o")n, (2-2) 

m—1 m =2 

where (-, -)n and (-, •) j m x n are the usual L 2 space and space-time inner-products, (■, -)/ m x n is the duality product 
between L 2 (/ m ; Tf _1 (n)) and L 2 (I m ; Hq (fi)). We note, that the first sum vanishes for u £ X JP. The dG(g) 
semidiscrete (in time) approximation u k £ X% of (TO) is defined as 

B(u k , <p k ) = (/, <Pk)ixn + (u 0 , for all (p k £ X%. (2.3) 

Rearranging the terms in (12.2b , we obtain an equivalent (dual) expression of B: 

M M—l 

= - '^2( u ’dt<p)imX n + (Vu, Vip) Ix n - Y ( M m. Mm)n + n- (2.4) 

m—l m =1 

Next we dehne the fully discrete approximation. For h £ (0. /v. 0 ]; ho > 0, let T denote a quasi-uniform 
triangulation of fl with mesh size h, i.e., T = {r} is a partition of f 1 into cells (triangles or tetrahedrons) r of 
diameter h T such that for h = max T h T , 

diam(r) < h < C'|r|» r j Vr £ T. 

Let Vh be the set of all functions in H q) (O) that are polynomials of degree r £ N on each r, i.e. V k is the usual 
space of conforming finite elements. To obtain the fully discrete approximation we consider the space-time finite 
element space 

x k,h = {vkh £ L 2 (I;V h ) \v kh \i m £V q {V h ), m = 1,2,..., M } , q> 0, r > 1. (2.5) 

We define a fully discrete cG(r)dG(q) solution u kk £ X?T by 

B(u k fn (p k h ) {f , ^Pkh)I xfi T (uq, ^P kk )ft for <tll ^p k h £ ^ k .h • 


(2.6) 
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2.1. Main results. Now we state our main results. 

2.1.1. Global pointwise best approximation error estimates. The first result shows best approximation 
property of cG(r)dG(q) Galerkin solution in L°°{I x O.j norm. For N = 2 and q = 0, r = 1, the result can be 
found in lf28l for convex polygonal domains. A similar result showing optimal error estimate is established in (9j, 
Thm. 1.2. We are not aware of any pointwise best approximation type results for N = 3. 

Theorem 2.1 (Global best approximation). Let u and Ukh satisfy O) and (lT6t respectively. Then, there 
exists a constant C independent of k and h such that 

T 

||u - u k h\\L°°(ixn) < Gin— |ln/i| inf ||u - xlU°°(/xfi)- 

K xexf h 

The proof of this theorem is given in Section[6] 

2.1.2. Interior pointwise best approximation error estimates. For the error at the point xo we can obtain 
a sharper result, that shows more localized behavior of the error at a fixed point. For elliptic problems similar 
results were obtained in Il34l 1361. We denote by Bd = Bd{x o) the ball of radius d centered at Xq. 

Theorem 2.2 (Interior best approximation). Let u and Ukh satisfy (O and m, respectively and let 
d > 4 h. Let t £ I m with some m £ {1,2 ,M} and Bd CC Cl, then there exists a constant C independent of 
h, k, and d such that 

T \ 

\{u-u kh ){t,x 0 )\ < Gin — |ln h\ inf l \\u - xlU~((o,t m )xB d (x 0 )) 

K X&X k ’ h 

+ ^ll u - xlU“((0,t m );L 2 (n)) + h\\v(u - x)l|L-((o,t m ) ; i 2 (n)) 

The proof of this theorem is given in Section|7] 

3. Elliptic estimates in weighted norms. In this section we collect some estimates for the finite element 
discretization of elliptic problems in weighted norms on convex polyhedral domains mainly taken from E). 
These results will be used in the following sections within the proofs of Theorem l2. H and Theorem l2.2l 

Let xo £ O be a fixed (but arbitrary) point. Associated with this point we introduce a smoothed Delta function 
ll36l Appendix], which we will denote by S = 5 Xo . This function is supported in one cell, which is denoted by t X() 
and satisfies 



(xJ)t X0 = X(zo), V X £V r (T X0 ). (3.1) 

In addition we also have 

||«|| w ..p(n) < Ch— N ^\ l<p<oo, s = 0,1. (3.2) 

Thus in particular H^Hl 1 ^) < G, ||^||z. 2 (f 2 ) < Ch ~and ||^||l=o(o) < Ch~ N . Next we introduce a weight 
function 

a(x) = i/|x - x 0 | 2 + K 2 h 2 , (3.3) 

where I\ > 0 is a sufficiently large constant. One can easily check that a satisfies the following properties: 

lk~^llL 2 (n) < G|lnft|5, (3.4a) 

|Vcr| < C, (3.4b) 

|V 2 (j| < G|(j —1 (3.4c) 

max cr < G min cr, Vr. (3.4d) 

For the finite element space 14 we will utilize the if projection Ph : L 2 (Ll) —f 14 defined by 


(PhV, x)n = ( v, x)n, Vx e V h 


(3.5) 
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the Ritz projection Tf ,: Hq (f2) —> V), defined by 

(Vi?^,V X )n = (Vt;,Vx)n, V X 6 V h , (3.6) 

and the usual nodal interpolation ih : Co(fl) —>■ Vh- Moreover we introduce the discrete Laplace operator 
A h : Vh —>• Vh defined by 

{~&hVh,x)n = Vx)n, Vx £ V/,. (3.7) 

The following lemma is a superapproximation result in weighted norms. 

Lemma 3.1 (Lemma 2.3 in ITTl ). Let Vh £ Vh- Then the following estimates hold for any a, ft £ R and K 
large enough: 

\\a a (ld-i h ){a i:! v h )\\ L 2 (n) + h\\o a W{ld -ih)^ v h )\\ L 2 {n) < c/i||er“ + ' 3 ~ 1 tv l || i 2 (n) , (3.8) 

\\a a (Id-P h )((J 0 v h )\\ LH n) +h\\a a V(Id-P h )(aPv h )\\ L 2 ( n) < ch\\o a+ f i - i v h \\ L , ( n ) . (3.9) 


The next lemma describes a connection between the regularized Delta functional 5 and the weight a. 

Lemma 3.2. There holds 

N — N ~ N ~ 

llc^lk^n) + /i||aTV<S|| L 2 (n) + 11 o'11 i, 2 (n) < C. (3.10) 

The proof of the above lemma for N = 2, for example, can be found in [|9) and for N = 3 in (77), Lemma 2.4. 
The next result shows that the Ritz projection is almost stable in L°° norm. 

LEMMA 3.3. There exists a constant C > 0 independent on h, such that for any v £ L°°(n) n 

||-R?i^||L“(n) < C , |ln/i|||u|| i =o( 0 ). 

For smooth domains such result was established in 031 . for polygonal domains in Bill , and for convex polyhedral 
domains in G3 Thm. 3.1]. In the case of smooth domains or for convex polygonal domains the logarithmic factor 
can be removed for higher than piecewise linear order elements, i.e. r > 2. The question of log-free stability 
result for convex polyhedral domains is still open. 

Next lemma is rather peculiar and can be thought as weighted Gagliardo-Nirenberg interpolation inequality. 
The proof is in El , Lemma 2.5. 

Lemma 3.4. Let N = 3. There exists a constant C independent of K and h such that for any f £ Hq(Q), 
any a, 0 £ R with a > —^ and any 1 < p < oo, ~ + y = 1 holds: 

\W a f\\ 2 mn) < C , ||tr“-^/|| i p(n)|| t 7“+ 1 +^V/|| LP , (n) , 

provided ||cr“ _/3 /|| iP (Q) and ||o' Q+1+/3 V/|| i px n ) are bounded. 

4. Weighted resolvent estimates. In this section we will prove weighted resolvent estimates in two and 
three space dimensions. We will require such estimates to derive smoothing type estimates in the weighted norms 
in Section 0 Since in this section (only) we will be dealing with complex valued function spaces, we need to 
modify the definition of the L 2 -inner product as 

(■ u,v)n= / u(x)v{x)dx, 

Jn 

where v is the complex conjugate of v and the finite element space as V/, = V), + iV) t . 

In the continuous case for Lipschitz domains the following result was shown in [39]: For any 7 £ (0. there 
exists a constant C independent of z such that 

\\{z + A)~ 1 u|| lp( o) < 1 l|tt|| L p ( n), z£C\E 7 , 1 < p < 00 , v € L p ( f2), (4.1) 

where S 7 is defined by 


E 7 = { z £ C |arg z\ < 7 } . 


(4.2) 
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In the finite element setting, it is also known that for any 7 G (0, f) there exists a constant C independent of h 
and 2 such that 

||(z +Aft)~ 1 x||z/»(fi) < Y^-pjllxllL-(n), for 2 G C \ E 7 , forallx'GV?,. (4.3) 

For smooth domains such result is established in [O and for convex polyhedral domains with a constant containing 
|ln h | in fl7l . In |20l the above resolvent result is established for convex polyhedral domains for some 7 G (0, 77 ), 
but with a constant C independent of h. 

Our goal in this section is to establish the following resolvent estimate in the weighted norm. 

THEOREM 4.1. For any 7 G (0, -|), there exists a constant C independent ofh and z such that 

||cr ^(2 +A ft ) _ 1 x|| L 2 (n) < C 1 1 n /t l || 0 -f x || Z|2(n); forz G C\E 7 , 

for all x G V/,, where E 7 is defined in (14.2I >. 

4.1. Proof of Theorem 14. II for N — 2. For an arbitrary 7 G ¥/, we define 

u h = (z + A h )~ 1 X, 

or equivalently 

z(u h , ip) - ( Vu h , V+) = (x, ip), V+ G Y h . (4.4) 

In this section the norm || • || will stand for || • ||i 2 (n). To estimate 11 ouh || we consider the expression 

HcrVu/jll 2 = (V(ct 2 u? 1 ), Vu h ) - 2(oVau h , \7u h ). (4.5) 

By taking ip = — Ph(a 2 Uh ) in (14.4b and adding it to (14.5k we obtain 

- z\\au h \\ 2 + \\aVu h \\ 2 = F, (4.6) 

where 

F = +1 + F 2 + F 3 := ~(cr 2 u h ,x) + (V(cr 2 u? l - P h {cr 2 u h )),X7u h ) - 2 (crS7ou h ,Vu h ). 

Since 7 < |arg 2 | < n, this equation is of the form 

e la a + b=f, with a,b> 0 , 0 < |a| < 7 r — 7 , 

by multiplying it by and taking real parts, we have 

a + 6 <(cos(|)) I/I < (sin (D) \f\ = C y \f\. 

From (14.61 ) we therefore conclude 

MII <TU k|| 2 + HcrVztft.il 2 < Cj\F\, for 2 G C\E 7 . 

Using the Cauchy-Schwarz inequality and the arithmetic-geometric mean inequality we obtain, 

1*11 = |(cr 2 Ztft, x)| < l|crztft,|| ||cr X || < CC 1 \z\^ || C 7 X || 2 + ^-\\(TU h \\ 2 . 

To estimate F 2 we use Lemma f+Tl the Cauchy-Schwarz and the arithmetic-geometric mean inequalities, 

1*21 < ||cr —1 V(cr 2 ztft, — P/ i (fT 2 u/ l ))||||( tVu/i|| < -^||aV^|| 2 +CC 7 ||^|| 2 . 

Finally, using the properties of 0 , we obtain 

|* 3 | < C\\u h \\\\aVu h \\ < -±-\\ a Vu h \\ 2 + CCJu h \\ 2 . 
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Combining estimates for F[s and kicking back, we obtain 

Mlk^ll 2 + ||<tVu/i|| 2 < c 2 {\z\- l \\*x\\ 2 + IKI| 2 ) • (4-7) 

Thus, in order to establish the desired weighted resolvent estimate, we need to show 

K|| 2 <C|ln/ l | 2 M- 1 || ( 7 X || 2 . (4.8) 

To accomplish that, testing (14.4b with ip = Uh, we obtain similarly as above 

|z||K|| 2 + ||Vu4 2 <C 7 |/|, for ^€<C\E 7 , 
where / = ( Xi u h )■ Using the discrete Sobolev inequality (see lf33l Lemma 1.1]), 

\\vh\\L°°(Q) < C\\nh\ 2 \\Vv h \\ L 2 {n) , \/vh £ 14 , 
and using the property of a (I3.4ab . we obtain 

klll«h|| 2 + l|Vu/J 2 < ^11(7x111,2^)||o-" 1 u ft || L2(n) 

< C 7 11 (7x142(0) 11 (7 1 1| r, 2 (n) ||tr/i||/, oo (Q) 

< C 7 |ln/i|||CTxl|z, 2 (n) II Vttft 142(0) 

< (7 2 |ln/i| 2 ||(7x|| 2 2(o) + ^||Vm? 1 || 2 2 (o)- 

Kicking back 41| Vm^II 2 2 (op we establish (14.81 ) and hence Theorem l4.1l in the case of N = 2. 

4.2. Proof of Theorem 14. II for N — 3. The three dimensional case is more involved and we require some 
auxiliary results. For a given point xq € U, we introduce the adjoint regularized Green’s function G = G x ° (x, z ) 
by 

G = G Xo (x,z) = (z + A)~ 1 5 
and its discrete analog Gh = G^°(x, z) £ Yh by 

G h = Gl°(x,z) = (z + A h )- 1 P h d, 

which we can write in the weak form as 

z('p,G h )-(V<p,VG h ) = {<p,$), Vp>£Y h . (4.9) 

From da we have the following result. 

Lemma 4.2 (El). Let Gh £ Yh be defined by ( 14.9b . There holds 

||G fc |4a(o) <C\\nh\i. 


Lemma 4.3. Let Wh £ Yh be the solution of 

z(w h ,<p) - (Vw,,, Y<p) = (/, ip), Vtp£ Y h 
for some f £ Li (LI). There exists a constant C > 0 such that 

IKI4~(fi) < c 'l lnft l^l/ll L i(o)- 


Proof. There holds 


w h (x 0 ) = z(wh, Gh) - ( Yw h ,VG h ) = ( f,G h )• 
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Hence, 

\w h (x 0 )\ = \(f,G h )\<\\f\\ Li{u) \\G h \\ L3{n) . 

Applying Lemma l4~2l we obtain the result. □ 

LEMMA 4.4. Let Vh G Yh be the solution of 

z(v h , <p) - (Vv h , Vip) = (/, <p), ftp G Y h , 

and f G L 1 (H). There exists a constant C > 0 such that 

IKIU 3 (Q) < C'|ln/i| 3 ||/|Ui(n)- 


Proof. We consider a dual solution Wh G ¥/, defined by 

z(<p, w h ) - (Vip, Yw h ) = (ip, v h \v h \), V<p G Y h . 

There holds 

IWIi 3 (n) = z(v h ,w h ) - ( Vv h ,Vw h ) = (f,w h ) < ||/||Li(n)ll^lli”(n)- 

By Lemma l4~3l that also holds for the adjoint problem, we have 

IKIU°°(fi) < C\^h\i\\v h \v h \\\ L 3 in) < C\lnh\^\\v h \\ 2 L3{U y 

Thus, we get 

lkft|li 3 (n) — C|l n ^l 3 ll/l|L 1 (f 2 )ll u /i|li, 3 (n)- 

Canceling ll'f/ill^Q) completes the proof. □ 

With these results we proceed with the proof of Theorem l4.1l for N = 3. 

Proof. For an arbitrary \ G Yh we define 

u h = (z + A h^x- 

or equivalently 

z{uh,<p) - (Yu h ,Y<p) = (x,<p), V<p€Y h . 

To estimate ||cr 5 u^|| we consider the expression 

\\a^Yu h \\ 2 = (V (a 3 u h ), Yu h ) - 3 {a 2 Yau h , Yu h ). 

By taking p = — Ph(a 3 Uh ) in (14.1 Ob and adding to 114.1 11 1 , we obtain 

- z\\<J%U h \\ 2 + ||<T^ Vltft || 2 = F, 

where 

F = Fx +F 2 + F 3 := -( P h {a 3 u h ),x ) + {Y(a 3 u h - P h (a 3 u h )),Yu h ) - 3(a 2 Yau h ,Yu h ). 
Since 7 < |arg z\ < 7 r, this equation is of the form 

e la a + b=f, with a,b> 0 , 0 < |a| < 7 r — 7 , 

by multiplying it by e~ ‘r and taking real parts, we have 

a + 6 <(cos(!)) I/I < (sin(|)) \f\=C 7 \f\. 


(4.10) 

(4.11) 

(4.12) 
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From (14.12b we therefore conclude 

| 0 |||ai^|| 2 + ||a«V^|| 2 <C' 7 |F| ! for 2 e<C\£ 7 . 

Using the Cauchy-Schwarz inequality and the arithmetic-geometric mean inequality we obtain, 

|Fi| = \(a 3 u h ,x)\ < Muhllllo^xll < C'C' 7 |2;|" 1 ||a5xl| 2 + t^rlk^ll 2 - 

To estimate F 2 we use Lemma lXTl the Cauchy-Schwarz and the arithmetic-geometric mean inequalities, 
\F 2 \ < \\a-?V(a 3 u h -P h (a 3 u h ))\\\\aiVu h \\ < Vu^ 2 + CC^unf. 

Finally, using the properties of < 7 , we obtain 

|F 3 | < C\\aiu h \\\\a^u h \\ < -i-||of Vu/i || 2 + CC^\\a^u h \\ 2 . 

Combining the estimates for F- s and kicking back, we obtain 

\z\\W^u h \\ 2 + llcrfvztft.ll 2 < C (N^lMxll 2 + W^unf) . 

Thus, in order to establish the desired weighted resolvent estimate, we need to show 

lltTiztftH 2 < C'|ln/t| 2 |^|- 1 ||crfx|| 2 . 

To accomplish that, we consider the expression 

—^||crf ztft. || 2 + ||crf Vztft|| 2 = —z(u h , au h ) + (Vu h ,X7(au h )) - (Vu hl 'Vau h ). 

Testing (14.1 Ob with <p = Ph(cruh ) we obtain similarly as above 

|^|||crfztft,|| 2 + ||cTfVztft,|| 2 < Cyl/I, for zeC\S 7 , 

where 

/ = fi + h + h ’■= -{Ph(vu h ),x) + (V(ctu/, - P h (au h )), Vu h ) - (V cru h ,Vu h ). 

Using the Cauchy-Schwarz inequality and the arithmetic-geometric mean inequality, we obtain 

l/ll = |(crztft.,x)| < 11 cr f ztft, 1111 cr f X11 < ^ || cr“ f ztft. || 2 + ^||crf x|| 2 - 
To estimate f 2 we use Lemma [3~TI the Cauchy-Schwarz and the arithmetic-geometric mean inequalities, 
I /2 1 < ||cr _ ^ V(crzt/ 1 — /^(crztft.)) || ||crf Vztft || < ■^-\\(T^Wu h \\ 2 + CC 1 \\(T~^u h \\ 2 . 
Finally, using the properties of cr, we obtain 

I/ 3 I < C'llcr - f Ztft, || ||crf Vztft. || < -4_ ||crf Vztft || 2 + CC 1 \\(J~^U h \\ 2 . 

Combining estimates for f- .s and kicking back, we obtain 

\z\\W^ u hf + ||crf Vztft,|| 2 < C ^ 11 cr f ztft 11 2 + ||crf xll 2 ) • 

To estimate ||cr — f ztft,|| we use Lemma l3~4l with a = /3 = —h and p = 3, to obtain 

||0—f Ztft.|| < C'||M/ l ||/ 3(n) ||VM ft ||%^. 


(4.13) 

(4.14) 


(4.15) 


(4.16) 
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Using Lemma l4~4l we have 

IMIl3 ( q) < C|ln/i|3||x||z,i(o) < C'|lnft.|3||(7 _ 5|||| (T l x || < C'|ln/i|s|| t T§ x ||. 

To estimate II we proceed by the Holder inequality 

ll v «h|| £ i (n) <C\\nh\*\\aiVu h \\ L2{n) . (4.17) 

Thus, using (14. 15b and the above estimates, we have 

Ml W^Uh\\ 2 + ||cr^Vzt/i|| 2 < C (||« h |U» ( n)||V« A || i3(n) + ||crf x || 2 ) 

< C (|ln/t|||c 7 ^Vu h ||||cr^xll + IM^xll 2 ) 

< C|ln/i| 2 ||cr^x|| 2 + Vu/jII 2 . 

Kicking back ||er'JVit/ l || 2 , we finally obtain 

lk^,|| 2 <C'|ln/ l | 2 N- 1 ||<7i X || 2 , 
which shows (14. 1 4b and hence the theorem. □ 

5. Maximal parabolic and smoothing estimates. In this section we state some smoothing and stability 
results for homogeneous and inhomogeneous problems that are central in establishing our main results. Since we 
apply the following results for different norms on 14, namely, for IJ'(fA) and weighted L 2 (H) norms, we state 
them for a general norm |||-|||. 

Let |||-||| be a norm on 14 (extended in a straightforward way to a norm on Vh) such that for some 7 £ (0, f) 
the following resolvent estimate holds, 

|||0 + A /l ) _ 1 x||| < jylllxlll, for z £ C \ £ 7 , (5.1) 

M 

for all x G V/,, where £ 7 is defined in (14.21 and the constant Mh is independent of z. 

This assumption is fulfilled for |||-||| = |H|iP(n), 1 < p < 00 , with a constant Mh < C independent of h , 
see 12JJ| , and for |||-||| = ||<7^ • 11z. 2 (£ 2 ) with Mh < C|ln/i|, see Theorem l4.ll 

5.1. Smoothing estimates for the homogeneous problem in Banach spaces. First, we consider the homo¬ 
geneous heat equation (fO-i.e. with / = 0 and its discrete approximation Ukh G Tff’jT defined by 

B{v,khiTkh) Tkh,0^ ^Pkh G -^k,h' (5.2) 

The first result is a smoothing type estimate, see M Theorem 13], cf. also M Thmeorem 5.1] for the case of 
the L 2 norm. 

Lemma 5.1 (Fully discrete homogeneous smoothing estimate). Let |||-||| be a norm on Vh fulfilling the 
resolvent estimate (ED- Let Ukh be the solution of ED. Then, there exists a constant C independent of k and h 
such that 

sup |||0 t Ufe/ l (f)||| + sup |||A h u fcfc (t)||| + Av; 1 i|[«fc/ l ] rri _i||| < ^111^4ttolll, 

te/m telm t rn 


for m = 1,2,..., M. For m — 1 the jump term is understood as [ufch]o = u^ h 0 — J4uo- 

For the proofs of Theorem 12.11 and Theorem 12.21 we will need an additional stability result, which is also 
formulated for a general norm |||-||| fulfilling (15.1k 

LEMMA 5.2. Let |||-||| be a norm on Vh fulfilling the resolvent estimate ED- Let Ukh be the solution of ED- 
Then there exists a constant C independent ofk and h such that 

M 

m—1 



\\\d t Ukh(t)\\\dt - 


\\\A h Ukh(t)\\\dt + \\\[u kh ] r 


< CM h In I P h u 0 HI. 
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For m = 1 the jump term is understood as [ukh ]o = u^h o — PhUo- 
Proof Using the above smoothing result, we have 


M 

E 


\\\dtu kh (t)\\\dt + / \\\A h u k h{t)\\\dt + |||['U^]m-i||| 


M 


m= 1 


— E ( Slip (t) HI + SUp HI AjiUkhit) III H" krri III [^/ch]m— 


tei , 


telrr 


M 


< CM h Y, T^lWPhUolW < CM h lnj\\\P h u 0 \l 


where in the last step we used that Y^m=i ^ ^ T- ^ 

5.2. Discrete maximal parabolic estimates for the inhomogeneous problem in Banach spaces. Now, we 
consider the inhomogeneous heat equation dm. with uq = 0 and its discrete approximation Ukh € defined 

by 


B {ttkh, T kh) — (f)Tkh'), ^Tkh G ^k,h‘ (5-3) 

The following discrete maximal parabolic regularity result is taken from [ 16j Theorem 14], 

Lemma 5.3 (Discrete maximal parabolic regularity). Let |||-||| be a norm on Vh fulfilling the resolvent esti¬ 
mate <EB and let 1 < s < oo. Let Ukh be a solution of dm. Then , there exists a constant C independent of k 
and h such that 



\\\dtu k h{t)\\\ s dt 



\\\A h u k h(t)\\\ s dt 


M 


+ E [Ukh] m —11 


\m= 1 


T 


< CM h In — [ / \\\P h f(t)\\\ s dt) , 


with obvious notation change in the case of s = oo. For m = 1 the jump term is understood as = u^h o- 

REMARK 5.4. As mentioned above the assumption dm is fulfilled for |||-||| = ||-||LP(n) an d any 1 < p < oo 
with Mh < C and for |||-||| = ||tJ^ • 11 i 2 (f 2 ) with Mh < C|ln/i|. Therefore the results of Lemma \5.1\ Lemma 13721 
and Lemma\5f3\are fulfilled for these two choices of norms with the corresponding constants M 

6 . Proof of Theorem 12.11 Let t £ (0, T] and let Xo £ V. be an arbitrary but fixed point. Without loss of 
generality we assume t £ (tM-i, T\. We consider two cases: t = T and tM-i < t < T. 

Case l,t — T: To establish our result we will estimate Ukh (T, xf) by using a duality argument. First, we 
define g to be a solution to the following backward parabolic problem 


—d t g(t,x) — A g(t,x) =0 (t,x) £ I x LI, 

g{t,x) = 0, ( t,x)£lxdLl , (6.1) 

g{T,x) = S x o, x £ Ll, 


where 5 = S Xo is the smoothed Dirac function introduced in dm Let gkh G be the corresponding 

cG(r)dG(q) solution defined by 


B(pkh,gkh) = Pkh(T,x o) Vp k h£X%’' h . (6.2) 

Then using that cG(r)dG(g) method is consistent, we have 

u k h{T,x 0 ) = B{u k h,gkh) = B(u,g k h) 

M M—l 

= ^ ^ dtQkh)ImXQ H - (VU, S7gkh)lxO, ^ ^ (^m? [9kh] m)f2 H - (6* 3) 

m—l m—l 

= J\ + J '2 + J'.i + 1/4. 
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Using the Holder inequality we have 


M 


Ji < E IMU“Umxn)||difffc/t||Li(/ m ;Li(n)) 

m= 1 

M 

< IMU“(Ixn) lli 1 (J'm;i 1 (f2))■ 


(6.4) 


m= 1 


For J 2 we obtain using the stability of the Ritz projection in L°°(H) norm on polygonal and polyhedral domains, 
see Lem ma [Q1 


J 2 = (VRhU,X7gkh)ixn = ~{RhU, A h g k h)ixn 

< \\Rhu\\ L oc( Ix ty || A h g kh \\ L i(1-^(0.)) 

< Cjlnft,|Hullioo^xQ)|| A h g kh \\L^(i-L^(Q.)) 


(6.5) 


For ./.j and J 4 we obtain 


M -1 


M—l 


J3 A 5 ' ll u m||z/ 00 (n) II [*7fc/t]m ||l 1 (0) — ll w llL oc, (7xn) ^ ' II ||L 1 (n) I 

m=l m—l 

J 4 < ||u(T , )||i~(n)||fl , fc/ li MlUi(n) ^ H u llL=°(/xn)ll5fe /l)M ||Li(n)- 


( 6 . 6 ) 


Combining the estimates for J\, J 2 , J 3 , and J 4 and applying Lemma [Ql with |||-||| = || • ||li(o) and < C, cf. 
Remark BT41 we have 


M 


\m= 1 


\ukh(T,xo)\ < C'|lnft|||u|| Il oQ(/ X Q) ( 5^ ||^s , fcfe||i i (/m;L 1 (n)) + WAhgkhW^ii-^in)) 

m- 1 n 

+ E \\[gkh]m\W{Q) + IK-^mIU 1 ^) 


< C'|lnft|ln^||u|| i oo ( j xn) ||P h 5|| I ,i ( n) 

T 

< Cjlnft|ln — ||u|| L «.( /x a), 

where in the last step we used the stability of the L 2 projection l\ with respect to the I? (fl) norm, see, e. g., Q 
and the fact that ||< 5 ||z, 1 (f 2 ) < C. 

Using that the cG(r)dG(g) method is invariant on X ?C, by replacing u and u k h with u — x and Ukh — X f° r 
an Y X C Xk'h , a nd using the triangle inequality we obtain 

\u(T, x Q ) - u kh (T, x Q )\ <Cln-Enft| inf ||u-xlli~(/xf 2 )- 

Case 2, t M -i < t < T: 

In this case we consider the following regularized Green’s function 


-d t g(t,x) - Ag(t,x) = 6 X0 (x)9(t) (t,x) elxil, 

g(t,x) = 0 , (t, x) £ I x dfl, 

g(T,x) = 0, x £ n, 

where 9 £ C' 1 (I) is the regularized Delta function in time with properties 


(6.7) 


and 


supp(9 c (t M -i,T), ||0||li(j m ) < C 

(8,<Pk)l M = <Pk(t), V(/7fe £ Vq{lM )■ 
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Let gkh be cG(r)dG(g) approximation of g , i.e. 

B(ipkh->9 9kh ) = 0 Pkh £ ^k]h m 
Then, using that cG(r)dG(g) method is consistent, we have 


V'kh(t) «^o) — {^kh 5 ^xo— B^Ukhigkh') — B{u,CJkh) 

M M 

= — y: (u, d t gkh)i m xsi + (Vm, Vgfc/,)zxn — Tl K) [gfcfe]m)a i 

m=l m=l 

where in the sum with jumps we included the last term by setting <jkh,M+i = 0 and defining consequently 
[Sfch]M = —gkh,M- Similarly to the estimates of Ji, J2, J3 above, using the stability of the Ritz projection in L°° 
norm on polyhedral domains, see Lemma tOI we have 

M M 

Ukh{t,x 0 ) = - ^2(u,d t g k h)ixn + (Vu, Vg k h)ixn ~ ^{um, [gkh\m)n 

m =1 m =1 

( M M 

T. + ||^7iPfch||i 1 (J;i 1 (n)) + ll[5fe?i]m|U 1 (0) 

m—1 m=l 

Using the discrete maximal parabolic regularity result from Lemma [531 with |||-||| = || • ||x,i(fj) and Mh < C, cf. 
Remark [5.4l we obtain 

T T 

Ukh(t,x 0 ) < C'ln-|ln/i|||u|| I/TO ( /xn) ||P^4 0 || z ,i (n) ||6'|| L i (jM) < Cln —|ln/i|||u|| L oo (/xn) . 

As in the first case this implies 

T 

\u(t,x 0 ) - u k h(t,x 0 )\ < C\n—\lnh\ inf ||u - xlU“(/xn)- 

K xexl'h 

This completes the proof of the theorem. 

7. Proof of Theorem l2.2l To obtain the interior estimate we introduce a smooth cut-off function u> with the 
properties that 



w( x) = 1, x £ Bd 

(7.1a) 

u(x) =0, x G fl \ B 2 d 

(7.1b) 

Vw < Cd~\ |V 2 w| < Cd~ 2 , 

(7.1c) 


where Bd = Bd(x 0 ) is a ball of radius d centered at £ 0 . 

As in the proof of Theorem l2.1l we consider two cases: t = T and tM-i < t < T. In the first case we obtain 


u kh (T,x 0 ) = B(u kh ,gkh) = B(u, g k h) = B(oju,g kh ) + B(( 1 - u>)u,g kh ), (7.2) 

where g is the solution of (l6Tt and g k h G AT?T is the solution of ( 16.2l i. The first term can be estimated using 
the global result from Theorem l2.ll To this end we introduce u = um and the cG(r)dG(g) solution u kk G 
debited by 


B(u kh - u, ip kh ) = 0 for all ip kh G X^. 


There holds 

T T 

B(u, g kh ) = B(u kh ,g kh ) = u kh (T, x 0 ) < Gin — |ln/i|||u|| i <*>(j x n) < Gin — |ln/i|||it|| L =o( /xB2d ). 

This results in 


T 

\u k h(T,x 0 )\ < Gin— |ln/i|||u|| L =o( 7xB2<i ) + B(( 1 -u>)u,g kh ). 


(7.3) 
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It remains to estimate the term B((l— ui)u, g k h)- Using the dual expression (12.4b of the bilinear form B we obtain 

M 

B(( 1 -u})u,g k h) = - ^((1 - u) u ,d t gkh)i m xn + (V((l - uj)u), Vg kh )i x n 

m= 1 
M—1 

w )«rn, [gkh]m)n + ((1 - w)w(T), g kh>M )n 

m= 1 

M (7.4) 

E N N V / 

(o’ - " 2 "(1 - u})u,a^d t gkh)i m xu + (V((l - u)u),Vg k h)ixn 

m= 1 
M—1 

E at iv iv at _ 

(cr~^(l -w)u m ,cr^[ 5fcft ] m )n + (cr - ^(l - u)u{T), g~ h M )u 

m= 1 

= ■/] + j'2 + 2s + J4. 

For Ji, using that <r~^ < Cd~% on supp(l - u) C O \ Bd and (1 — ui) < 1, we obtain 

M 

Ji < ||a-T(l-a;)u||L- (/;L 2 ( n )) \W^d t g k h\\L^(i m -L^(Q.)) 

m=1 (7.5) 

M v ' 

N \ N 

< Cd 2 ||u||L“(/;L 2 (n)) 2_^ Ho -2 dtgkhW^ilm-^L 2 ^))- 

m =1 

To estimate J 2 , we define U (1 oj)u and proceed using the Ritz projection R k defined by (13.6b . There holds 

(Vip(t), Vg kh (t))n = (V R h ip(t),Vg k h(t))a. = -(R h ip(t),A h g kh {t))n 

= ~(R h ip{t), A h g k h(t)) Bd/2 - {Rhip(t), A h g k h(t))n\ Bd/2 
< \\Rh'4>(t)\\ L °°(B d/2 ) \\Ahgkh(t)\\ L i(B d/2 ) 

+ Cd~~ \\Rh4>(t)\\ L 2(n\ Bd/2 )\\(T— A h g k h{t)\\ L 2(n\B d/2 ) 

< \\Rh4’{t)\\L™{B d/2 )\\A h g k h(t)\\Li-(a.) + Cd~^\\Rh4>{t)\\ L 2^\\a^A h g k h(t)\\ L 2(Q,), 

where we used < Cd~% on U \ B d / 2 . In the interior pointwise error estimates | [36l Thm. 1.1] with F = 0, 
choosing x = 0, s = 0, g = 2 and using the triangle inequality and the fact that supp (/>((:) C f2 \ B d , we have 

\\R h m\\L~(B d/2 ) < C\lnh\\m)\\ L -(B d ) + Cd-^\\R h m\\L H a) = Cd-Z\\R h m\\L> { n). 

Using a standard elliptic estimate and recalling ip = (1 — uj)u we have 

\\Rhi>(t) |U 2 (ft) < II^WIU 2 (n) + II ip(t) ~ Rhi/j{t) |U 2 (fi) 

< IIV’WIU 2 ^) + c/i||vV'(f)||i2 (n) 

< ||M(i)IU 2 (n) + ch\ |(1 - w)Vu(t) - VccM(f)|| i 2 (n) 

< c||u(i)||i 2 (n) + c/i||Vu(f)|| L 2 (n) , 

where in the last step we used |Vw| < Cd -1 < C/i -1 . 

Therefore we obtain 

(Vip(t), V< 7 fc/j(i))n < Cd~^ (|H0IU 2 (n) + c/i||Vu(f)IU 2 (n)) + || f 7^A^ feft ,(f)|| i 2 (n) ^ . 

This results in 

Ji < Cd 2 (||u||L“(/;L 2 (n)) + cB||VM|| i oo( /;i 2 ( n ))) ^||AftPfcft||x,i(J ; £,i(fi)) + Ik 2 A/ l pfc/ l || i i(/ ; L 2 ( fl ))^ . 

(7.6) 
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For J3, similarly to Ji we obtain 


M—l 

N \ N 

J 3 < ||<7 2 (1 — w)w||L”(/;L 2 (n)) 2_^ Ik 2 [ 9kh\m\\L 2 (Q .) 

m—l 

M—l 

jv * ■y n 

< Cd 2 ||u|| L oo(j ;i 2( n )) 2^ Ik 2 kfc/t]m||i 2 (Q). 
m=l 


(7.7) 


Finally, 


Ji<Cd 2 |k|| L =o( 7 ;L 2 (n ))|k 2 g kh jM IU 2 (n)- 


(7.8) 


wc nave 


Combining the estimates for J\, J 2 , J 3 , and J 4 , 

B((l - a ))v,g kh ) < Cd~% (|k|U=°(/;L 2 (fi)) + ch\\ Vw|| L <» (/;L 2 (n)) ) 

/ M 

/ % ^ iv iV 

X I 2_^ lk~k3fcti||L 1 (/ m ;L 2 (n)) + || A/ l gfckUi(/;L 1 (Q)) + |k~ X h g k h\\ l^(I-L 2 (Q.)) 

M—l ' 

y lk T [ffkfcUMn) + lk 2 <?fcMrlU 2 (n) 

m=l 


n— 1 


+ 


For the term || XhgkhWL 1 ^^ 1 ^)) we apply Lemma 15721 with |||'||| = || ■ ||li(q) and Mh < C and for all weighted 
terms with |||-||| = ||<r t (-)||L 2 (fi) and Mh < Cjlnft.|, cf. Remark [L4l resulting in 

B(( 1 — u)v,gkh) < Cd In —|ln/i| (IML”!/;^ 2 ^)) + /i|l Vw||L“(/;L 2 (n))) (||-P/k|Ui(n) + |k^~ Ak||.L 2 (n)) 

N T 

<Cd~ In —|ln/i| (|M|L°°(/ ; L 2 (fi)) + (i||VM|| i oo(/ ;i 2 ( n ))) , 

where in the last step we again used the stability of the L 2 projection with respect to the L 1 norm, the fact that 
Ikllz, Rfi) < C, and Lem ma im for the term || cr 2 PhS\\^(n)- Inserting this inequality into (17.31) , we obtain 

\u k h{T,x 0 )\ < Cln—|ln/i| (jk|U“(Jx.B 2d ) + d~~ (||w||i,~(/;z, 2 (n)) + h\\ Vu|| L «>(/ ; L 2 (n)))) ■ 

Using that the cG(r)dG(q) method is invariant on Xf.' T h , by replacing u and Ukh with u — x and u kh — X f° r any 
x e XK, we obtain Theorem l2.2l for the case t = T. 

In the case (m-i < t < T we proceed as in the proof of Theorem [2T] using the dual problem (16.71) instead 
of O- Then, we proceed as in the above proof using in the last step the discrete maximal parabolic regularity 
from Lemma [531 instead of Lemma 15721 This completes the proof. 
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